function [ xdot ] = model(t, x)
%MODEL Summary of this function goes here
%   Detailed explanation goes here

global alpha beta reinforcementsFnc;

xdot = zeros(2, 1);
xdot(1) = -alpha * x(2) + reinforcementsFnc(t);
xdot(2) = -beta * x(1);

end

